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Abstract 

Simultaneous intensity measurements of protons, helium, and heavy nuclei 
can be used to describe the effects on the source energy spectrum of solar 
modulation and diffusion occurring in source regions and in interstellar space. 
The energy spectra of heavy nuclei (Z > 10) have been determined in the energy 
region 200 - 700 Mev/nucleon for the year 1961 and compared with similar data 
obtained for protons, a-particles, and heavy nuclei during the years 1961 and 
1963. These results have indicated that the observed intensities could be 
adequately explained by the following: (a) the modulating mechanism as de- 
scribed by a model given by Parker compared with the electric deceleration 
model of modulation, (b) a source differential rigidity or total energy spectrum 
rather than a kinetic energy spectrum having an exponent of 2.35 (as derived 
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from experimental data available in the high energy range 10 10 -10 14 eV/nucleon, 
that is little affected by solar modulation) and (c) a path length of interstellar 
matter traversed that is decreasing or constant with decreasing energy for 
energies £15 GeV/nucleon. 

Introduction 

The energy spectra of the low energy cosmic ray nuclei observed at the 
earth depend on (a) the energy spectrum of nuclei emanating from the source 
(b) the propagation of these nuclei in the interstellar space (diffusion and ac- 
celeration) and (c) the modulation of the energy spectrum in the vicinity of the 
solar system. These various effects could be delineated from a systematic 
study of the simultaneous measurement of the energy spectra and intensity of 
various components of the cosmic radiation as a function of time. Protons and 
helium nuclei have the same energy loss per nucleon for particles of the same 
velocity, but have charge to mass (Z/A) ratios that differ by a factor of two. 
Therefore these two species of nuclei can be used to separate the rigidity and 
velocity effects of solar modulation, a-particles and heavy nuclei have a Z/A ratio 
that is almost 0.5, but have different rates of energy loss ( Z 2 //3 2 c 2 where /3c is 
the velocity of the nucleus); thus their relative intensities are similarly affected 
by solar modulation, but depend to a large extent on their diffusion in interstellar 
space and source regions. 

Systematic studies of the differential energy and intensity spectra began in 
the year 1961 for the protons and helium nuclei components [see, e.g., Fichtel 
et al, 1964]; the earlier data on such energy spectra were obtained from de- 
tectors flown in balloons under a few g/cm 2 of atmospheric depth and are 
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subject to uncertainties in the extrapolation procedures adopted. The heavy 
nuclei component was recently studied by Lim and Fukui [l965] and Fichtel et 
al. [l965] using emulsions flown in satellites and rockets during the year 1963. 
We report here those observations made on the energy spectra of cosmic ray 
nuclei of charge Z >10 in the energy interval 200-700 MeV/n during the period 
1961. Measurements of the energy spectra for proton and helium nuclei for 
this period have been previously made by Fichtel et al. [l964]. This paper is 
referred to hereafter as paper I. These observations are combined with the 
present ones and with similar spectra obtained during 1963 for the heavy nuclei, 
a -particles, and protons in order to expand on the three aspects mentioned 
previously. 

Experimental Details 

A stack of nuclear emulsions exposed to a mean altitude of ^2.4 g/cm 2 of 
residual atmosphere in a balloon near Fort Churchill during 1961 for a period 
of ~10 hours is used for the present purpose. The details regarding the stack 
size, the flight curves, description, and orientation of the stacks before and 
during the flight are described in paper I. 

The G-5 emulsions in this stack were scanned along a line parallel to the 
top edge and 5 mm inside the stack. The scanning line had a length of about 
5 cm. Tracks having a projected zenith angle dp < 50°, dip angle £29°, and 
ionization £64 I 0 (the minimum of ionization), were accepted. These tracks 
were followed through the stack till they stopped, interacted or, left the stack. 

The charge and energy of these nuclei were determined by a combination of 
S-ray counting and range and scattering measurements. The general pro- 
cedures adopted were discussed in previous literature in detail, Aizu et al., 



5 


[i960]. Since, in the present work, the main emphasis is on the determination 
of the energy spectrum, not much attention is paid to resolving with care in- 
dividual neighboring charges. We describe here briefly the methods adopted in 
the three cases. 


(a) Stopping tracks — The charges of the tracks were identified by making 
8-ray counts at the point where they crossed the scan line. Two S-ray counting 
procedures were used: short S-ray counting (in which S-rays extending be- 
yond a projected distance of >3.2^ from the track on either side were counted). 
At the end of the range an integral S-ray count was performed over a distance 

of ^2000^ on favorable tracks that ended in G5 emulsions. These integral 
counts, for nuclei of different charges, are related by a "similarity law". [See 
Aizu et al., I960] For calibration purposes, tracks having known charge values 
in emulsions exposed in previous rocket flights [ Fichtel et al ., 1965] and having 
the same ionization characteristics were used. The energy of each nucleus was 
estimated from the range and charge of the nucleus. 

(b) Tracks leaving the stack or interacting — These were identified by 
counting short and long S-rays (at least 200 in number) at the scan line and at a 
point close to the point of exit (1 cm from the exit) or interaction point. The 


calibration plot was derived from the well-identified long stopping tracks de- 


scribed in (a) above, and using the relation, the 8 -ray density N s 



where a and b are constants and Z and /3c are charge and velocity of the nucleus. 


The energy was estimated either for the change in the 8-ray density with range 
or from scattering measurements. 


Each track was followed backwards to the edge of the stack and an energy 


spectrum of particles at the top of the stack was obtained. 
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The procedure adopted for extrapolating the observed spectrum to the top 
of the atmosphere was the same as that adopted by Aizu et al ., [i960]. The ef- 
fective mean vertical atmospheric depth is 2.4 g/cm 2 ; in view of the low value of 
the matter traversed, it is assumed here that all the fragmentation corresponding 
to 2.4 g/cm 2 of air occurs at the top of the atmosphere. Thus, the energy of 
each nucleus has been first corrected for ionization loss in air and then for loss 
due to fragmentation. The following parameters are used: 

^vh = 12.9 g/cm 2 , \ H = 16.5 g/cm 2 , P VHVH ~ 0.28 ± 0.28 

P vh ( sh) = 0.33 + 0.15, P SH( sh) = 0.15 ±0.05 (VH: Z> 20 and SH: Z = 10-19} 

Results 

Figure 1 shows the energy spectrum of the particles obtained at the top of 
the stack. The number of nuclei having a particular charge value is shown in 
the box corresponding to the energy interval. The differential intensities of 
H (Z >10) nuclei obtained for the top of the atmosphere as a function of energy 
are given in Table 1 and are plotted in Figure 2. Also shown in this figure are 
the differential energy spectra of proton and helium nuclei obtained in this stack 
and at the same time by Fichtel et al . [l964]. The Mt. Washington neutron 
monitor counting rate on the day of the flight is 2148. 

The energy spectra of P and a nuclei have been best studied for the year 
1963 from the low energy region ('■/L00 Me V/nucleon) to a fewGeV/n(Balasubrahmanyan 
and McDonald [l965], Fichtel et al. [l965B], Frier and Waddington [l965], Ormes 
and Webber [l964]). The spectra thus obtained are also shown in Figure 2. The low 
energy rocket data as obtained by Fichtel et al . [l964] are combined with the satellite 
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TABLE 1 

Differential intensities of heavy nuclei (Z> 10). 


Kinetic energy (MeV/n) 

Flux (dJ/dE) 

Particles/m 2 sec sv MeV. 

160 - 275 

0.0027 ± 0.0007 

275 - 400 

0.0035 ± 0.0008 

400 - 600 

0.0026 ± 0.0005 

600 - 700 

0.0019 ± 0.0006 


data of Lim and Fukui [l965] to obtain an integrated energy spectrum of H-nuclei 
in the low energy region. The energy spectrum obtained by Koshiba et al . [l963] 
during a period close to solar maximum in a balloon flight (altitude ~2 g/cm 2 ) is 
also shown in this figure. The neutron monitor data corresponding to the dates 
of the flight are indicated in the figure. It can be noted, from a comparison of 
the spectra as a function of time that (a) there is a decrease in intensity of nuclei 
as solar activity increases and (b) there is a shift of the maximum towards 
higher energy with increasing solar activity. There will be a further discussion 
of the experimental results in combination with the satellite and rocket data in 
the next section. 

The helium spectrum has been measured in the present stack and is 
given in paper I. Similar simultaneous intensity measurements of a-particles 
and heavy nuclei were made by various authors ( Aizu et al . [i960], Evans [l963], 
McDonald and Webber [l962], and Biswas [l965]). The values of T Ha computed 
from these measurements are given in Figure 3. It can be seen that the ratios 
do not change appreciably with time, thus indicating a charge independence of 
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the modulation process in the energy region considered for the nuclei of same 
Z/A as has been predicted by various theoretical models. 

Discussion 

Source Spectrum 

The integral flux values of cosmic ray nuclei determined for the highest 
possible known energy, from solar maximum to solar minimum, are near the 
equator, (Cutoff rigidity -16GV); within the limits of experimental uncertainties 
and errors, these values do not show 7 variation. Thus it is reasonable to assume 
that the flux values obtained for this energy represent those outside the solar 
system. The integral flux values of protons have been measured at energies 
-10 12 eV, (See Webber [l964 b]), these flux values are shown in Figure 4. 

The galactic cosmic ray spectra are normally represented as a power law 
either in total energy, W, the integral intensity, J w (>W) of nuclei being given 
by C w W _r , or in kinetic energy, E, the integral intensity J e (>E) being given by 
C e E~ y and in rigidity R, the integral intensity J r ( >R) being given by C r R -y . The 
total energy, W, the kinetic energy, E, are expressed in GeV/n and rigidity R 
in GV, and C w , C e , and C r are constants. The power law in total energy is pre- 
dicated by various acceleration mechanisms. The statistical acceleration 
process proposed by Fermi [l954], the betatron mechanism of acceleration 
[ Ginzburg and Syrovatsky , 196l], and the acceleration in hydromagnetic shock 
waves produced during supernova explosions [ Colgate and Johnson , I960] all 
lead to such an energy spectrum. Ginzburg and Syrovatsky [l96l] have recently 
shown that, on the basis of general arguments concerning the equi-partition of 
energy in a turbulent magnetized plasma of cosmic dimensions, a power law in 
total energy of cosmic rays is a fundamental property of such a plasma. 
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A power law in kinetic energy has been used by Gloeckler , [l965] and Bal- 
asubrahmanyan et al. [l965] to explain the measured intensities of protons and 
helium nuclei by balloons and the IMP-I satellite; however, in their analyses 
they had to use a ratio of intensities of protons to helium nuclei (P/a) = 5.7, 
which is much lower than the value obtained near the equator, where there is 
very little modulation. For this reason, a source spectrum which is a power 
law in kinetic energy is not considered in the present analysis. 

The power law rigidity spectrum has recently been assumed by Kaplon 
and Skadron [l964] . There is, perhaps, no physical justification for assuming 
such a spectrum except for the fact that the solar flare spectrum could be 
represented by a rigidity spectrum. We assume in our analysis that the spec- 
trum in the source regiohs is either (a) a power law in total energy, (herein 
referred to as source total energy spectrum) or (b) a power law in rigidity, R 
(herein referred to as source rigidity spectrum). We further assume that all the 
components, protons, a-particles and heavy nuclei have the same exponent, y. 
This exponent is given as 1.35 from Figure 4. The values of constants C w and 
C f obtained for various nuclei are given in this figure. These values given 
here for protons may be compared with the values of C p = 5450 ± 80 and 
y p = 1.43 ± 0.10 (or 1.31 as deduced by Waddington , [i960]) obtained by Mc- 
Donald [l 9 5 8] in the energy range 4-15 GeV/n from a comparison of the flux 
values measured at Texas and Guam. 

Propagatio n in the Interstellar Medium 

Various models, e.g., the regular model and the diffusion model, [ Ginzburg 
and Syrovatsky, 196l] have been postulated for the diffusion of cosmic rays in 
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interstellar space and a model of confinement of these nuclei in source regions 
has been given by Kaplon and Skadron [l964] ; these models predict different 
variations of the path length with energy. We assume here the regular model of 
diffusion of cosmic rays in which all particles of a unique energy move along a 
definite path, for example, along the magnetic lines of force. We further assume 
that no acceleration takes place in the interstellar medium; the existing experi- 
mental evidence seems to indicate that a Fermi statistical acceleration mecha- 
nism cannot operate effectively in the interstellar medium [See, for example, 
Durgaprasad, 1965]. 

In Traversing a path length, X g/cm 2 of interstellar space, two competing 
processes occur, namely, the absorption of nuclei of the same species and the 
feeding-in of nuclei from heavier group through nuclear collisions and the loss 
of energy of particles due to ionization. For heavy nuclei of charge Z ^10, the 
energy loss and absorption of nuclei have been taken into account. Helium 
nuclei can arise from spallation of medium and heavy nuclei, besides being ab- 
sorbed; however, the intensities of these nuclei are so low that to a first approxi- 
mation, one could neglect the contribution from these nuclei. Consequently, the 
differential intensities of heavy nuclei (Z > 10) and alpha particles (Z > 2), 
(neglecting the contribution from the fragmentation of heavy nuclei and medium 
nuclei to alpha particles), can then be written as: 



and 
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Where (dJ/dR) ; and (dJ/dW). are the intensities of differential rigidity and 
total energy of i-nuclei, R. and W. are the rigidity and total energy of i-nuclei, 
DRSSI, DWSSI are the rigidity and total energy intervals at the solar system for 
the corresponding intervals DRSI and DWSI at the source region, after traversal 
through Xg/cm 2 of interstellar medium (for computing ionization loss, the 
medium is supposed to consist of neutral hydrogen) and A ; is the absorption mean 
free path of the i -nucleus. Using the values of C. given in Figure 4 and the 
absorption mean free paths, AH = 6.0 g/cm 2 , Aa = 36.0 g/cm 2 and Ap = 72.0 
g/cm 2 , the unmodulated differential intensities of heavy nuclei and a-particles 
for energy and rigidity spectra at the vicinity of the solar system have been 
calculated for various values of X; from these values the ratios of intensities, 
r Ha , were computed and are shown in Figure 3. The errors shown are derived 
from the errors on C.. These are compared with the available experimental 
data on ^ Ha in the energy region, 200-700 MeV/n, (Figure 3). It can be seen from 
this comparison that a matter traversal of ^2 g/cm 2 for the source rigidity 
spectrum as well as for the source total energy spectrum is consistent with all 
the experimental data. This result is also in good agreement with the recent 
experimental evidence on the observed ratios of He 3 /(He 3 +He 4 ) [Hildebrand et al ., 
1964] and of light L(Z = 3-5) to S-nuclei (Z> 6) [Hildebrand and Silberberg, 

1964]. However, it is not in agreement with the observations that suggest an 
increase in the matter traversed with decreasing energy [e.g., Apparao, 1961 
and Webber , 1965]. 

Propa gation in the Solar System 

The spectra thus obtained are modulated to a considerable extent by solar 
activity; different mechanisms have been postulated by means of which such 
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changes can occur and, of these, two models seem to be currently in vogue; 
Parker's model and the electric field model. Recently Fichtel, Durgaprasad and 
Guss [l965B] , from a study of the differential intensities of protons and 
a-particles, have shown that, assuming a similar spectra at the source, the gen- 
eralized Parker's model can adequately represent the intensity variations during 
a solar cycle; however, the source spectra had peculiar shapes and could not be 
represented by a unique power law. Frier and Waddington [l965], on the other 
hand, assumed a total energy spectrum with an exponent of 1.5 and tried to show 
that an electric field model could be used to predict these intensity variations. 

We consider here both these models in great detail for the two kinds of spectra 
assumed, and discuss the significance of these parameters. 

Electric field model: In this model it is assumed that the earth is at a 
positive potential and that positive particles do lose energy in traversal through 
the solar system; also they occupy a different volume in phase space, the density 
D in phase space being constant. Consequently, if P, M, and W are the momentum, 
mass, and total energy of these particles and J(P) and J(W) the corresponding 
differential intensities, the following relations hold good: 


= J(P)M = J(W) 

p3 P 2 


Constant 


( 2 ) 


Let us designate these intensities as being outside the solar system. Then the 
modulated intensities g(P)' and g(W)' could be obtained from: 


g(P')M = J (P)M 

( P .) 3 P 3 


g(W')M = J(W)M 
(W')3 W 3 


and 


( 3 ) 
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Where the momentum P' and total energy W' of the particles after the modula- 
tion is related to the quantities P and W before modulation from the relation 


W' 


W 


zev 



(4) 


Ze being the charge and V the positive electric potential of the earth with respect 
to infinity. Combining (1) and (3), we have the unmodulated intensity g(W) re- 
lated to the source spectrum from the relation: 


dj 

g(W) = dE 


HI - fa sill 

("1) ! - (m„c 2 ) 2 


C 


(», * a) 


2.35 



DWSI 

DWSSI 


) 


(5) 


[See, for example, Frier and Waddington , 1965] in the case of total energy 
spectra and 


dj 

g( R ) ~ dR “ 


W 2 - (m 0 c 2 ) 


2\2 


[_(W + A ) 2 - (m 0 c 2 )' 


3/2 


c. 


-x/Ai 


R 


2.35 


DRS I 
DRSSI 


( 6 ) 


in the case of total rigidity spectra. Values of g(W)and g(R) have been calcu- 
lated for different values of V and are compared to the flux values of protons 
and helium nuclei obtained for the year 1963 for a traversal of path lengths 
varying from 1 to 3 g/cm 2 . These are given in Figure 5. From this figure it 
can be seen that a good fit to proton data of the year 1963 can be obtained using 
a value of V~300 MV, but the agreement with helium nuclei data is not quite as 
good; this departure becomes more marked for the year 1961; a value of V = 500 
MV can be fitted to the proton data but the measured intensity values of helium 
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nuclei of the same potential for 1961 fall well below the calculated curve. In 
the case of source rigidity spectra, the proton data could be fitted using a value 
of V = 1000 MV whereas the helium nuclei data fall well below the line. Thus, 
besides the question concerning the existence of such high potentials of the order 
of a few hundreds of MV in the solar system, the experimental data could not 
be explained using this mode of modulation and the energy spectrum given above. 
We will now examine the applicability of the generalized Parker model to the 
present problem. 

Parker's model : In this model, the differential intensity of particles in the 
vicinity of the earth, g(i ), is related to the unmodulated intensity, J(i ), in the 
vicinity of the solar system as [Parker, 1963] 


Bi(i) 



(?) 


where i refers to the type of nucleus under consideration, V is the solar wind 
velocity, /3c the particle velocity, r is the distance from the sun and X is the 
scattering mean free path. By assuming that V/X is independent of r from r e 
to r 0 and zero thereafter, (7) can be written as 


g(i) 


J ( i ) exp 



where D(R) - 


3V ( r 0 ~ r e) 

MR) 


For X, Parker [l956] suggested the form 


X = 



Constant 
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which is equivalent to 


1_ _ Constant 

K (l + r 2 /r 0 2 ) 


where p is the radius of curvature of the particle of rigidity R, L is the diameter 

2BL 

of the scattering center, R 0 — and Bis the magnetic field strength. Thus, 

(7) can be rewritten as 


g(i) 


J ( i ) exp 


r k//3 (l + R 2 /R 0 2 )] 


( 8 ) 


k and R 0 being constants that depend on the solar wind velocity and the dimen- 
sions of the scattering centers. 

Thus the observed intensity at the earth is related to the source spectrum 
from the relation, 




exp 


-k//3 (l +R 2 /R 0 2 ) 



( 9 ) 


in the case of source rigidity spectrum, and from the relation 



S / DWSI \ 
W 2 . 35 VdwssiJ 


exp - k//3 ( 1 + R 2 / R q 2 ) e“ x/Ai 


(10) 


in the case of source total energy spectrum. For particles of a unique energy, 
there are three unknowns, namely K, R 0 , andX; and the experimental intensities 
have been measured for three nuclei, namely protons, helium, and heavy nuclei. 
The measured intensities of (dJ/dR) 0 and (dJ/dE) 0 for the years 1963 and 1961 
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are given in Figure 6 (a) and (b) respectively and these fluxes have been fitted 
for various values of K and R 0 as shown in these figures. 

It can be seen that a reasonable fit to the experimental points can be obtained 
in the case of source rigidity spectrum for values of k = 2.1 and R 0 = 1.6 for 
the year 1963 and for values of k = 2.2 andR 0 = 2.0 for the year 1961. A 
similar attempt has been made in the case of total energy spectra; the agreement 
is not so good as can be seen from the figure. The values of k and R 0 obtained 
here are of the order of magnitude that one would expect from the known solar 
wind velocities and expected sizes of magnetic irregularities [ Parker , 1963]. 

Thus one can conclude that the modulation mechanism described by Parker's 
model and a source rigidity (or a total energy) spectrum would adequately pre- 
dict the intensity variations with solar activity; the matter traversed at low 
energies being ~2 g/cm 2 . This value should be compared to a value of 
2.5 g/cm 2 of hydrogen used for nuclei of relativistic energies; thus, the value 
is either constant or else it decreases with decreasing energy. 

A slight decrease in the matter traversed by the radiation with decreasing 
energy is not quite unexpected, for if the low energy and high energy particles 
are produced at the same time, have the same age t (this is true if the bulk of 
cosmic ray nuclei of all energies are produced at the same time in the past, 
e.g., in a violent explosion of the nucleus of our galaxy or some other), and 
traverse regions of similar density (n), then from the simple relation, X = /3c nt 
where /3c is the velocity of the particle, one could infer that low energy particles 
traverse less matter than the high energy ones. This assumes a simple dif- 
fusion process and the energy or the rigidity spectrum that is observed in cosmic 
rays could be mainly dominated by such a diffusion process in the source regions 



17 


and in interstellar space. The solar cosmic ray events do show a rigidity de- 
pendent spectrum [e.g., Webber, 1964] and it is reasonable to expect a similar 
situation to occur in the galactic or extragalactic cosmic rays as well. Further 
studies of the energy spectrum and the intensity of the low energy cosmic ray 
particles using space vehicles would likely reveal some interesting information 
on the diffusion of particles in space and on the source spectrum of these 
nuclei. 
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Figure Captions 

Figure 1: Energy spectrum of heavy nuclei (Z >10) observed at the top of the 
stack. The numbers indicated in the box refer to the measured 
charge of the nucleus. 

Figure 2: Differential intensities of protons, helium, and heavy nuclei measured 
in 1963, 1961 and 1959 as a function of kinetic energy (MeV/nucleon) 
and rigidity. The numbers in brackets refer to the counting rate of 
the Mt. Washington neutron monitor. 

Figure 3: Ratios of intensities of heavy nuclei to alpha particles plotted as a 
function of Mt. Washington neutron monitor rate. Also shown in the 
text are the ratios predicted from Equation (1) (See text). The data 
shown refer to: A: Aizu et al . [i960]. P: present work. M: McDonald 
and Webber [l962], E. Evans [l963] and B. Biswas [l965]. 

Figure 4: Integral intensities of protons plotted as a function of kinetic energy 
(and rigidity) for E ~ 15 GeV where the effects of solar modulation are 
assumed to be negligible. The constants C p , C a and C H are those that 
define the integral intensities at the source (see text). The data shown 
are derived from Webber’s review article ( Webber , 1964B). 
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Figure 5(a): Differential spectrums of protons and helium nuclei measured in 
1963 plotted as a function of magnetic rigidity. The curves show 
predicted spectra for modulating potentials of 900 and 1000 MV, 
and for different values of interstellar matter traversed (given in 
brackets) . 

Figure 5(b): Differential intensities of protons and helium nuclei measured in 
1961 and 1963. The curves show predicted spectra for modu- 
lating potentials of 200, 300, 500 and 700 MV (see text: electric 
field model) and for different values of interstellar matter traversed 
(given in brackets). 

Figure 6 : Differential spectra of protons , helium and heavy nuclei measured 

in 1963 and 1961 plotted as a function of (a) magnetic rigidity (MV) 
and (b) kinetic energy MeV/n. The curves shown refer to various 
values of K and R 0 (see text: Parker's model) and X; the inter- 
stellar matter traversed as indicated in the figure. 
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